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TECHNICAL MEMORANDUM X-53961 


HOLOGRAPHIC OPTICAL DATA PROCESSING 
SUMMARY 


To process optical data is to perform mathematical operations, such as 
addition and multiplication on optical inputs. Using holography, these opera- 
tions may be performed on optical inputs having complex amplitudes. This 
report presents the theory and techniques necessary for the application of 
holography to optical subtraction and pattern recognition. 

The parameters affecting the successful recording of a hologram are 
discussed in detail. These include the photographic process, the coherence of 
the light source, the vibrational stability of the optical system, and the resolv- 
ing power of the recording medium. 

In optical subtraction it is shown how, by double exposure and the intro- 
duction of a 180-degree phase delay, a wave front may be subtracted from 
another. 

In pattern recognition, it is shown how a matched filter may be synthe- 
sized by holographic methods. Graphs of the sensitivity of the filtered output 
power to rotation of the input transparency are given. The appendixes show 
mathematically that a lens displays the Fourier spectrum of an input in its 
front focal plane on its back focal plane and how the filtered output is input 
translation invariant. 


INTRODUCTION TO HOLOGRAPHY 
Definition 


Holography is a method by which a highly complex wave front can be 
recorded and later reconstructed w|.th identical phase and amplitude from the 
recording medium. The recording is accomplished by photographing the inter- 
ference pattern created by a simple wave front, normally plane or spherical, 



interfering with the wave front to be recorded (Fig. 1). The recorded inter- 
ference pattern is unique in that only a wave front of identical phase and ampli- 
tude could interfere with the same simple wave front to recreate this pattern. 
This simple wave front is usually called the "reference wave front"; the wave 
front to be recorded is spoken of as the "signal wave front. " 


BEAM 



Figure i. Configuration for recording the wave front 
reflected from a three-dimensional object. 

Reconstruction of the recorded wave front is accomplished by reillumin- 
ation of the photographic plate with the same reference wave front used in the 
recording process. The photographic plate behaves as a diffraction grating, 
diffracting part of the reference wave front into a reconstruction of the stored 
wavefront (Fig. 2). 
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Figure 2. Configuration for reconstructing the wave front 
reflected from the three-dimensional object. 


The Photographic Process 


Exposure, development, and object beam to reference beam intensity 
ratio must be such that the response of the emulsion, measured in amplitude 
transmittance, is linearly proportional to the intensity being recorded. Unfor- 
tunately, emulsions do not have an entirely linear response to energy. Shown 
in Figure 3 is a typical transmittance versus exposure curve where amplitude 
transmittance is defined as the ratio of the amplitude of the light incident upon 
the plate to the amplitude of the light transmitted by the plate; the exposure is 
in units of energy so that it is proportional to intensity. If the emulsion is to 
have a linear response, the intensity modulation must be confined to the linear 
section of this curve. 
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Let the complex amplitude distribution (at the plate) of the reference 
wave front and the signal wave front be represented by.R(x,y) and S(x,y) 
respectively. If these wave fronts interfere at the surface of a square-law 
detector, such as a photographic plate, an intensity distribution, I(x,y), will 
be recorded, where 

I (x,y) = |R(x,y) +S(x,y)! 2 = |R(x,y)| 2 + |S(x,y)| 2 
+ R* (x,y)S(x,y) + R(x,y)S* (x,y). 

If the real amplitudes of the reference and signal wave fronts are approximately 
constant, we may regard the first two terms of the expression as constant. 

Then, the last two terms are the only modulation terms. 

Accordingly, the amplitude transmittance, T , of the plate will also 

3 , 

have a constant part, T , and a fluctuating part, f(x,y). Thus, 

T' (x,y) =T + f(x,y). 

a O 

To stay on the linear section of the curve, f (x,y) must be much less than T . 

This is satisfied if we pick |S(x,y)| much less than |R (x,y)| so that the inten- 
sity modulation terms are much less than the constant terms: 

R* (x,y)S(x,y) + R(x,y)S* (x,y) « |R(x,y)| 2 + |S(x,y)| 2 

In addition, if we pick R(x,y) such that 

|jR(x,y)| 2 + |S(x,y)| 2 J =*|R(x,y)l = 

where (E , T ) is a point that falls in the middle of the linear section of the 
o o 

curve (E being the desired exposure) , then transmittance modulation function 
f (x,y) will be linearly proportional to the intensity modulation. Therefore, 


E o V /2 


T a (x,y) = T q + )8t [r* (x,y)S(x,y) + R(x,y)S* (x,y)] , 

where j3 is the slope of the curve at (E , T ) and t is the exposure time. 

o o 
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If the plate is reilluminated with the same reference wave front used in 
the recording, the light transmitted will be 

R (x,y)T (x,y) = R (x,y)T + /3t|R(x,y)| 2 S(x,y) + /3tR 2 (x,y)S* (x,y). 

3. O 

The second term, /3t|R(x,y) | 2 S(x,y), is a reconstruction of the signal 
wavefront except for a simple constant (assuming the reference wave front has 
a constant amplitude distribution in the x,y plate) multiplied times its ampli- 
tude. This term creates the virtual image shown in Figure 2. The first and 
third terms do not contribute to the reconstruction of the signal wave front 
since they travel in a different direction. (Note that if we use R* (x,y) 
instead of R(x,y) to reconstruct the hologram, the third term becomes 
/3t|R(x,y) | 2 S*(x,y), which creates a real image. ) 


Coherence 


Since no source of light emits an infinitely long wave train, even highly 
monochromatic laser light has associated with it a temporal coherence length, 
AL, which is proportional to the duration. At, of a single wave train [1], 


AL = cAt — c/Av = 




AA 




where Av is the effective frequency range of the Fourier spectrum and A q is 

the mean wavelength. This distance determines the maximum path difference 
between two wave fronts from a common source that will give satisfactory 
fringe contrast. If the path length difference between the reference and the 
signal wave fronts exceeds the coherence length of the laser, a hologram will 
not be recorded. In practice, the coherence length must be experimentally 
determined for each light source since laser cavity expansion because of heat- 
ing will ultimately limit the theoretical value. 

It is also important to have a spatially coherent reference wave front. 
This can be accomplished by focusing the reference beam through a pinhole 
before it strikes the plate. The diameter of the pinhole should be on the order 
of the size of detail desired to be resolved in the image [2], It is important 
to note that the spot size of the beam at its focal point should be less than or 
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equal to the diameter of the pinhole if the system is to be diffraction limited, 
d 2 = SS = i. 22 Af/D, 

where d is the diameter of the pinhole, SS is the spot size, A is the wavelength, 
f is the focal length of the focusing lens, and D is the effective diameter of the 
lens. 


Stability 


The relative phase of the two interfering beams must be held constant 
during the exposure if the interference pattern is to be recorded. Hence the 
path length difference between the reference wave front and the signal wave 
front must be held constant to within a fraction of a wavelength of the light used 
to record the hologram. This is usually accomplished by either rigidly mount- 
ing all components of the optical system to a common heavy support which is 
isolated from floor vibrations by air mounts (or some cushioning material) or 
by using a high-powered pulse ruby laser to shorten the necessary exposure 
time to nanoseconds. 


Emulsion Resolving Power 


If the spacing between the interference fringes becomes smaller than the 
grain size of the emulsion, the fringes will not be recorded. Since the angle 
between the signal and reference wave fronts determines the spacing, this angle 
must not exceed a certain critical angle , 9 , determined by 

0 = 2 arc sin (A/2p), 

where A is the wavelength of the recording light and p is the resolving power of 
the emulsion. 
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APPLICATIONS 


Optical Subtraction 


Holography can be used to optically subtract one complex wave front 
from another. If the wave fronts are identical except for an additional term in 
one, the additional term can be emphasized by subtracting the two wave fronts. 

The technique is to superimpose two holograms on the same photo- 
graphic plate by double exposure. If the same reference wave front is used in 
each exposure, the reconstruction of the signal wave front will be the sum of the 
signal wave fronts used in each exposure. To make the wave fronts subtract 
rather than add, a 180-degree phase delay is introduced into the signal wave 
front to be subtracted. 

If S and S + AS are the complex amplitude distributions of the signal 
wave fronts for the first and second exposures, respectively, and R is the 
complex amplitude distribution of the reference wave fronts for both exposures 
(omitting the coordinate notation for abbreviation) , the light transmitted by the 
doubly exposed hologram will be 

RT = RT + fit [R 2 S + R 2 S*] + fi't' [R 2 (S + AS) + R 2 (S + AS)*] . 
a o 

(For multiple exposure holography, it is desirable to prebias the emulsion to 
a transmittance T^ at the beginning of the linear section of the T-E curve. Then 

the reference-to-signal amplitude ratio can be unity and the successive exposures 
can be chosen such that fi^^ i 3 2 t 2 = /3 3 t 3 ) 

If a 180-degree phase delay is introduced into the wave front to be sub- 
tracted and if fit ~ fi 1 1' , we can write 

RT = RT + fit f|R | 2 S exp (in) + R 2 S* exp (nr) + R 2 (S+AS) 
a o L 

+ R 2 (S+AS)*1 =RT +/3t|R| 2 AS + /3tR 2 AS*. 

-I o 

Effectively, S has been erased, leaving only a reconstruction of AS. 
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This was first done with Fourier transform holography [3] (Fig. 4). 
Because the Fourier spectrum of a signal input created by backlighting of a 
transparency is relatively insensitive to translation of the transparency (see 
Appendix A) , such a configuration is much easier to use than a system whose 
signal wave front is the backscatter off an object. However, in principle, it 
should be possible to subtract backscattered signals also. Figure 5 is a photo- 
graph of the holographic image resulting from the subtraction of the signal 
diffusely reflected by the right half of a small card from the signal diffusely 
reflected from the entire card (the right half of the card is much dimmer than 
the left half). The holographic configuration that was used. for this subtraction 
is shown in Figure 6. 



TRANSPARENCY WITH INPUT 
SIGNAL TRANSMITTANCE 


Figure 4. Configuration for optical subtraction by Fourier transform holography. 
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Pattern Recognition 


Pattern recognition can be 
accomplished by optical spatial filter- 
ing using holograms as matched filters. 
The purpose of such a system is to de- 
tect the presence of a given pattern by 
enhancing its optical signal and sup- 
pressing any extraneous signals caused 
by other patterns. This is accomplished 
by causing the desired signal to have a 
higher energy density at the surface of 
a detector than any of the extraneous 
signals. 


Figure 5. Photograph of the image 
reconstructed from a subtraction 
hologram made by the configuration 
shown in Figure 6 (The background 
for the number 69 was erased by 
optical subtraction.) 


The basic system for accom- 
plishing this task is shown in Figure 7. 
The pattern to be recognized is placed 
at Pi in the form of a transparency 
(such as a photographic negative). A 
plane parallel beam of coherent light is 


incident upon the transparency at Pj and 

is transmitted with complex amplitude s(x 1 ,y 1 ). The lens L 2 displays a complex 

amplitude distribution Sfc^Af, y^Af), at P 2 which is the Fourier transform of 

s (x ls yj). (See Appendix B: Fourier Transformation by A Lens. ) Placing a 

transparency at P 2 with complex amplitude transmittance T (x^Af , y^Af) 

a 

causes S(x 2 /Xf, y 2 /Af) to be multiplied by T (x 2 /Af, y^Af). Lens L 3 then dis- 

a 

plays the Fourier transform of this product at plane P 3 . 


If the signal s (x,y) is to be recognized, the transparency at P 2 should 

have 


T (x 2 Xf, y 2 Af) = kS* (x/Af, y 2 /Af) 
a 

so that 

T (x 2 /Af, y2 /Af)S(x 2 /Af, y 2 /Af) = k| ! S(x 2 /Af, y 2 /Af)l 2 . 
a 
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LASER 



ii 


Figure 6. Configuration for optical subtraction of backscattered signals. 




The curvature of the wave front incident at P 2 is thereby canceled, leaving a 
plane wave which can be focused to an approximate point on P 3 by lens L 3 as in 
Figure 7. Such a transparency is called a matched filter for input signal s (x,y). 

However, until the invention of holography, it was not possible to make 
a transparency with complex amplitude transmittance kS* (x 2 /Af, y 2 /Af) . With 
holography, a matched filter for any given input signal can be easily made. 

Shown in Figure 8 is the system for making the matched filter. To make the 
matched filter for an input signal, s (x,y) , s (x,y) is Fourier transformed onto 
plane P 2 , Also incident on plane P 2 is a plane parallel beam that has complex 
amplitude distribution R exp (i27rQ!yj) on P 2 , where a = sin 0/k and R = constant. 
If a photographic plate is placed at P 2 , a hologram can be recorded with com- 
plex amplitude transmittance 

T (x 2 /Af, y^Af) = R 2 + S 2 + RS exp (-i27TQ!y2) + RS^ exp (i27ra:y2). 

3 . 


The fourth term is the desired transmittance except for the linear phase term 
exp (i27ray 2 ). However, this term is desirable since it separates the desired 
transmittance from the other terms by inclining the desired transmittance at 
an angle, 6, off axis (Figure 9). Multiplying T (x 2 /Af, y^Af) by 

cl 

S(x 2 /Af, y 2 /Af) gives a reconstruction otf the plane parallel wave front incident 
upon P 2 except for the | S(x 2 Af, y^Af)! 2 multiplier. Then, lens L 3 will take 
the Fourier transform of this term which approximates a Dirac delta function 
centered at coordinates (O, aAf) on P 8 at which a light-sensitive detector can 
be located (Appendix B). 

To experimentally verify the theory, a matched filter was made for a 
pattern of vertical slits by the configuration shown in Figure 10 on a 4-by 5-inch 
Kodak 649-F photographic plate. To make the filter fairly selective, a stop 
one-half inch in diameter was located at the center of signals Fourier spectrum 
to act as a high-pass filter. This eliminated the dc term and the lower fre-r 
quency (spatial) diffraction. It was found that, as was expected mathemat- 
ically, the effect of translation of the input transparency in the (xj, y t ) plane 
was only an equal shift of the output point in the (-x 3 , -y 3 ) output plane. How- 
ever, the output was sensitive to rotation of the input transparency. Graphs of 
the sensitivity of the matched filter to rotation of the input transparency are 
plotted in Figures 11 and 12. 
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Figure 8. Configuration for making a matched filter by holography. 



Figure 9. Schematic showing off-axis displacement 
of a holographically constructed matched filters output. 
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OUTPUT 



Figure 10. Experimental pattern recognition configuration 
used to obtain results plotted in Figures i i and i 2. 




Figure ii. Graph of the response of the matched filters output 
(measured in units of power) to a rotation of the input 
transparency about its vertical axis. 
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Figure 12. Graph of the response of the matched filters output 
(measured in units of power) to a rotation of the input 
transparency about its longitudinal axis. 
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APPENDIX A 


SHIFT INVARIANCE OF MATCHED FILTER OUTPUT 

To find the effect on the filtered output of a shift of the input, let the 
input signal be shifted by an amount, a, such that s(x 1 ) becomes slxj-a). 
Then, the Fourier transform of sCxj-a) will be 


F s(xj-a) = f s(x 1 -a)exp (-2irif x^dxj. 

“* _oo 

Let Xj-a = x', then x t = x* + a and dx 4 = dx* so that, 


F^slxj-a) J = J s(x’) exp [(-27rif x ) (x' + a)]dx* 

oo 

= exp (- 27rif a) f s(x') exp (-2nif x 1 ) dx’ * exp (— 27rif a) F [s(x 1 )] 
x J x x ‘ 

_ OO 

= exp (-27rif a) S (f ) . 

xx 

Thus, the only effect in the frequency plane is a linear phase shift. 

Multiplying by the amplitude transmittance of the matched filter, 

T (f ) [exp ( - 27rif a) S(f )] =R|S(f )l 2 exp (-2?rif a) exp ( 2 ^ 0 x 2 ). p 
ax xx x x 
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To find the distribution on plane P 3s we must take the Fourier transform of 
this product: 


F [T (f ) exp (-27rif a)S(f )] 
ax x x 


OO 

= f R ! S (f )| 2 exp [(-27Tif ) (a-aAf)] exp (-27rif x 3 )df 

^ X X XX 

— oo 


where y 9 = f Af. If R |S(f ) I 2 =* K, which is a constant, then 

j & x x 


OO 

F [T (f ) exp (— 27rif a)S(f )] = K f exp [ (-2irif ) (x 3 +a-o:Af)]df 

3, X XX ^ X X 

_oo 


= K<5 (x 3 + a - a\f) , 

where <5(x 3 +a-aAf ) is a Dirac delta function centered at coordinate (-a+oAf) 
on the x 3 axis of plane P 3 . 

Thus the only effect of a shift in the (x^y*) plane is an equal shift in 
the (-x 3 , -y 3 ) plane. 
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APPENDIX B 


FOURIER TRANSFORMATION BY A LENS 1 


A thin lens acts as a phase transformation on an incident wave front. 
Referring to Figure B-l, the total phase delay may be written 

$(x,y) = kn A(x,y) + k[A 0 - A(x,y)] , 



Figure B-l. Schematic of a wave 
front with complex amplitude 

incident upon a lens and the complex 
amplitude UJ transmitted by the lens. 


where n is the index of refraction of the 
lens material and the index of the sur- 
rounding medium is unity. Hence, the 
transmitted wave front is related to the 
incident wave front by 

Uj (x,y) * exp ^j$(x,y)^U (x,y). 

The thickness function, A(x,y) , 
may be written 

A(x,y) = Aj(x,y) + A 2 (x,y) , 

and, by reference to Figure B-2, it is 
seen that 


At(x,y) = A^ - (Rj - ^Ri - x 2 - y 2 j = A^ - R t (l - 


A 2 (x,y) = A 0 ^ - (-R 2 - \/r! - x 2 - y 2 ^ = A 0 ^+ R 2 (i - 


1. This development follows the approach made by J. W. Goodman in his book, 

Introduction to Fourier Optics. 
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Figure B-2. Schematic for calculating the thickness 
function of a lens . 

If only paraxial rays are considered, then the following approximations 
can be made: 



Substitution of these approximations gives: 
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With the substitution of this expression, the lens phase transformation term can 
be written 

exp (j4(x,y)^ = exp (jknA 0 ) exp ~j“ (x 2 + y 2 ) , 

where f, the focal length, is defined as 




To find the distribution U_(x„,yJ in the back focal plane of the lens, the 

III 

Fresnel diffraction formula can be applied: 


U f (x f ,y f ) 


jAf 


// u i (x i l 3 'i ) exp i (x f - x i ) 


+ (y f - Yj) 2 j dXj dy J# 

Expanding the quadratic terms is the exponent 


U f (x f ,y f ) 


exp(.jkf) 


jm exp hi ( x t +y 0] 

i ~ (xf + yf)] exp [ (XjXj y^) ] dXj dyj ■ 
The expression for U| in terms of is 


exp 


UJ = exp ( j* (x^yf) Uj (x 1 ,y 1 ) = exp (jknA 0 ) exp [ -j^“ (x 2 + y 2 ) 
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Substitution causes the quadratic phase terms within the integral to vanish 
leaving 


U (x v ) = exeiMl ex n 

VW jkf exp 


j il( x f +y 0 


oo 

i ! /u i (x 1 ,y l ) e xp[-jf ( Vf+yi y ( ) 


dXjdy^ 


exp (jkf) 
jkf 


exp 


2z 


( x l + y|)] F [hW 


where the symbol F [ ] denotes the Fourier transform operation and the trans- 
form is evaluated at frequencies 



Consider a plane object with 

amplitude transmittance t (x ,y ) as 
^ o o o 

in Figure B-3. If a monochromatic 
plane wave of amplitude A is incident 
on the object, the disturbance trans- 
mitted will be 


U (x. ,y ) - A t (x ,y ). 
060 000 


Again, the Fresnel diffraction formu- 
la can be used, this time to relate 

U. to U : 

1 o 



Figure B-3. Configuration for 
Fourier transformation of a plane 
object amplitude transmittance by 
a lens. 


U i (x i-V = 


exp (jkf) 
jAf 


If u Jx >y J ex P 

J J 000 




2 


+ 



dx dy . 
o o 
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The form of the formula is that of a convolution of U (x ,y ) with 

o o o 

h(x,-x ,y -x ), where 
1 o 1 o 


exp (jkf) 

1 V J 1 J o' jAf 


h(x,-x_ ,y,-y_) = “ — exp \j 


_k 

2f 


<x r x o )!+(y i- y o ,! ] 
Using the convolution theorem, we can write: 


U^x^y^ 


F [WV] p [ h(x l' y 1 ) } 


Recalling our expression for U^.(x^,y^) and substituting for 

V x rV " 


where 


U f (x £ , y f) exp j-j_k (x> + y>)]p [u 0 (x 0 ,y 0 )] P [hlx^ 


F 


h(x i ,y i ) ] =zex P ( j kf) ex P [-fjr < x f + yp ]• 

h(Xj,y^) J and dropping the constant phase delay 
U ( (x ( ,y ( )=F[U o <X o ,y o )]. 


Substituting for F 
terms , it is found that 


It is therefore seen that a lens Fourier transforms an optical input in 
its front focal plane onto its back focal plane. 
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Mr. W. A. Darbro 
Mr. J. W. Parker 
Mr. G. Hintze 


S&E-SSL-PM 
Mr. E. Shriver 
Mr. P. Espy 
Mr. E. Klingman 
Mr. K. S. Clifton 
Mr. J. McGuire 
Mr. D. Jex 

S&E-SSL-PO 

Dr. S. G. Frary 
Dr. T. R. Edwards 
Mr. J. R. Williams 
Mr. W. W. Moore 
Mr. P. W. Tashbar 
Mr. R. W. Jenkins (20) 
Mr. R. L. Kurtz 
Dr. N. M. Reid 
Mr. J. Simmons 
Mr. N. Norden 
Mr. R. Freeman 

S&E-SSL-S 

Dr. W. H. Sieber 
Dr. O. K. Hudson 

S&E-SSL-ST 

Mr. B. S. Blake 

S&E-SSL-SG 

Mr. A. G. Loughead 

S&E-SSL-SE 

Mr. L. S. Yarbrough 
Mr. W. L Chisholm 

S&E-SSL-T 

Mr. W. Snoddy 
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INTERNAL (Concluded) 


EXTERNAL 


S&E-SSL-TE 

Mr. E. R. Miller 

S&E-SSL-TT 

Mr. B. P. Jones 
Mr. P. Craven 

S&E-SSL-TR 

Mr. G. M. Arnett 
Mr. J. Zwiener 

S&E-ASTR-R 

Mr. J. C. Taylor 
Mr. E. J. Reinbolt 
Mr. J. L. Randall 
Mr. R. B. Hoover 

R-QUAL-ARA 

Mr. M. C. Mcllwain 

R-QUAL-ARM 
Mr. J. B. Beal 
Mr. R. L. Brown 


A&TS-PAT 

Mr. L. D. Wofford, Jr. 
PM-PR-M 
A&TS-MS-IP (2) 
A&TS-MS-IL (8) 

A&TS-TU (6) 

A&TS-MS-H 

DIR 

DEP-T 


National Aeronautics & Space Admin- 
istration 

Washington, D. C. 20546 
Attn: Dr. H. H. Kurzweg RR 
Dr. W. A. Menzel REI 
Mr. J. W. Keller RV-i 
Mr. A Reetz RV-i 
Mr. D. Novik RFE 

Scientific and Technical Information 
Facility (25) 

P. O. Box 33 

College Park, Maryland 20740 

Attn: NASA Representative (S-AK-RKT) 

Mr. A. R. Shulman 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 

Mr. G. Bonner 
Bldg. 31 — Room 225 
Manned Spacecraft Center 
Houston, Texas 77058 

The University of Tennessee Space 
Institute 

Tullahoma, Tennessee 37388 
Attn: Dr. T. H. Gee (4) 

Dr. J. D. Trolinger 

Dr. M. K. Testerman 
University of Arkansas Graduate 
Institute of Technology 
P. O. Box 3017 
Little Rock, Arkansas 
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EXTERNAL (Concluded) 


Dr„ J e A. Muscari 
Martin Marietta Corporation 
P, O, Box 179 (Dept 1610) 
Denver, Colorado 

Mr. J. C. Pernicka 
Executive Electronics 
Stapleton International Airport 
Hanger #5 

Denver, Colorado 80207 



